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RETRIEVAL  OF  CIRRUS  RADIATIVE  AND  SPATIAL  PROPERTIES 
USING  INDEPENDENT  SATELLITE  DATA  ANALYSIS  TECHNIQUES 
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1.  INTRODUCTION 

The  role  of  clouds  in  climate  is  receiving  emphasis 
in  the  US  Global  Change  Research  Program.  The 
current  consensus  among  researchers  is  that  clouds 
have  a  net  cooling  effect  on  the  Earth's  climate,  but 
how  much  effect  is  still  in  question  (Arking,  1991). 
Important  factors  include  cloud  fraction,  altitude,  and 
optical  depth.  The  only  continuous  and  truly  global 
observations  of  clouds  are  provided  by  satellites; 
however,  questions  remain  concerning  what 
information  can  be  retrieved  and  how  representative  it 
is  of  true  cloud  conditions.  To  this  end,  the 
International  Satellite  Cloud  Climatology  Project 
(ISCCP)  is  a  worldwide  effort  to  obtain  accurate 
estimates  of  global  cloud  data  from  satellites  for 
climate  models  (Rossow  and  Schiffer,  1991). 

Cirrus  is  one  of  the  most  poorly  quantified  clouds. 
As  a  part  of  ISCCP,  intensive  observations  of  cirrus 
clouds  were  taken  in  the  autumn  of  1986  over 
Wisconsin.  During  this  First  ISCCP  Regional 
Experiment  Cirrus  Intensive  Field  Observation  (FIRE 
Cirrus  IFO),  coordinated  measurements  from  satellite, 
aircraft,  and  ground-based  platforms  were  made  of 
cirrus  clouds. 

This  paper  deals  with  the  verification  of  cirrus  cloud 
information,  both  spatial  and  radiative,  obtained  for  a 
1986  FIRE  cloud  scene  using  measurements  from  two 
independent  sensors  onboard  the  NOAA-9  polar 
orbiting  satellite.  In  addition  to  the  wide  variability  in 
properties  common  for  other  types  of  clouds,  cirrus 
clouds  have  the  added  complexity  of  transmissivity 
values  t that  span  the  entire  possible  domain  0  <  t  <  1 . 
Thus,  uncertainties  exist  in  thin  cirrus  cloud  amount, 
altitude,  thickness,  and  optical  properties  as  retrieved 
from  satellite  because  the  measured  cirrus  signal  is 
affected  additionally  by  an  unknown  radiation 
component  from  below. 


2.  CIRRUS  RETRIEVAL  MODEL  ATTRIBUTES 

Multiple  cirrus  analysis  models  are  useful  for 
verifying  independently  derived  cirrus  characteristics 
within  the  same  cloud  scene.  The  results  from  two 
such  models  are  compared  in  this  study;  the  C02 
Slicing  technique  and  a  multispectral  Advanced  Very 
High  Resolution  Radiometer  (AVHRR)  infrared  cirrus 
analysis  model  developed  by  d’Entremont  et  al.  (1990) 
called  the  AVHRR  Infra-Red  Cirrus  model  (AIRC). 
While  the  fundamental  requirement  to  detect  and 
analyze  thin  cirrus  is  the  same  for  both  models,  the 
capabilities  of  the  two  models  differ  in  numerous 
respects.  Table  1  lists  the  important  attributes  of  the 
HIRS  COz  and  AIRC  cirrus  retrieval  models.  The  most 
important  differences  are  the  spectral  bands  and 
spatial  resolutions  of  the  satellite-detected  radiances 
used  by  each  model. 

In  recent  years  the  C02  Slicing  technique  has  been 
extensively  applied  to  obtain  cloud  statistics  on  a 
global  scale  (Wylie,  Menzel,  and  Woolf,  1991).  The 
technique  estimates  cirrus  altitude  and  the  product  of 
cloud  fraction  (N)  and  emissivity  (e),  called  “effective 
emissivity  (Ne)”  using  radiances  from  the  NOAA 
High-Resolution  Infrared  Radiation  Sounder  (HIRS) 
C02  spectral  channels  in  the  13.4-14.2  pm  range.  It  is 
capable  of  detecting  transmissive  cirrus  clouds  at 
upper  tropospheric  altitudes  above  levels  where  the 
C02  channels’  weighting  functions  peak.  At  these 
HIRS  COz  wavelengths  the  differential  absorption  of 
atmospheric  radiation  allows  the  C02  Slicing  method 
to  detect  transmissive  clouds.  When  a  transmissive 
cloud  is  found,  the  HIRS  Channel  8  (11.1  pm)  IR 
window  brightness  temperature  measurement  Tg  is 
used  along  with  an  estimate  of  the  surface  skin 
temperature  Tsfc  to  compute  effective  emissivity  and 
cloud  altitude.  If  the  C02  channel  radiances  do  not 
detect  the  presence  of  transmissive  cloud,  then  an 
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Attribute 

AIRC 

C02  Slicing 

Input 

Satellite 

Data: 

AVHRR  LAC/GAC 
IR  Window  Chs  3, 
4,  and  5 

HIRS  CO2  Sounder 
Chs  4  -  7,  And 

IR  Window  Ch  8 

Analysis 

Resolution: 

Nominal: 

LAC  - 1  X  1  km2 
GAC -  4  X  4  km2 

Nominal: 

25  X  25  km2 

Retrieved 

Parameters: 

IR  “Effective” 
Emissivities; 
Effective  Altitude; 
Cirrus  Fraction 

IR  Effective 
Emissivity; 

Effective  Altitude 

Table  1.  Attributes  of  Cirrus  Analysis  Models 


emissivity  of  1.00  is  assumed  and  Tg  is  compared  to 
the  temperature  profile  to  assign  a  blackbody  cloud 
top  pressure.  Cloud  climatologies  produced  using  the 
C02  slicing  method  have  been  extensive,  focusing  on 
the  geographical,  seasonal,  and  diurnal  changes  of 
cloud  cover  (Wylie  and  Menzel,  1989;  Menzel,  Wylie, 
and  Strabala,  1991;  Wylie,  Menzel,  and  Woolf,  1991). 
Through  these  studies  the  first-order  radiative 
characteristics  of  transmissive  cirrus  clouds  are  now 
better  ascertained  on  a  global  scale. 

C02  Slicing  cloud  analyses  have  been  compared 
with  lidar  measurements  and  NWS  ground-based 
cloud  reports  (Wylie  and  Menzel,  1989).  Comparisons 
with  other  satellite  techniques,  however,  have  not  yet 
been  made.  This  paper  presents  the  early  results  of 
the  first  such  comparison  with  the  AIRC  cirrus  retrieval 
model.  AIRC  generates  cirrus  emissivities, 
transmissivities,  and  optical  depths,  along  with 
effective  cirrus  layer  altitude.  AIRC  determines  these 
cloud  properties  using  the  differences  among  AVHRR 
3.7,  10.7,  and  11.8  pm  IR  window  brightness 


Figure  la.  AVHRR  3.7  pm  Image  for  the  0930  UTC 
28  October  1986  FIRE  IFO 


temperature  measurements.  The  model  is  based  on 
simple  radiative  transfer  principles  that  account  for 
both  the  semi-transparent  nature  of  cirrus  clouds  and 
the  attenuative  effects  of  atmospheric  water  vapor  in 
the  AVHRR  thermal  window  regions. 

Resolution  differs  considerably  between  the  two 
models.  For  the  HIRS  data,  this  FOV  is  large: 
nominally  490  km2,  corresponding  to  a  25  km 
diameter.  It  is  expected  that  most  HIRS  FOVs  are  not 
uniformly  cloud  covered  and  contain  large  variations  in 
cloud  density.  This  has  the  effect  of  not  allowing  an 
analysis  of  the  HIRS-derived  cirrus  emissivity  itself;  it 
is  intrinsically  tied  to  cirrus  cloud  fraction.  Like  the 
C02  technique,  the  AIRC  model  can  only  resolve  the 
"effective”  emissivity.  However,  AIRC  is  applied  to 
Local  Area  Coverage  (LAC)  AVHRR  data  with  a 
nominal  FOV  of  0.8  km2,  corresponding  to  a  1  km 
diameter.  Compared  to  HIRS  pixels  the  AVHRR  LAC 
pixels  have  less  variation  in  cirrus  cloud  density  and, 
except  near  cirrus  edges,  are  more  often  either 
entirely  clear  or  entirely  cloudy.  This  dichotomous 
state  of  the  AIRC  cloud  amount  (N  =  0  or  1)  allows 
direct  determination  of  cirrus  emissivity. 

3.  DATA  ANALYSIS 

Nighttime  AVHRR  Local  Area  Coverage  (LAC)  and 
coincident  HIRS  C02  sounder  radiance  data  valid  at 
0930  UTC  on  28  October  1986  were  collocated  on  a 
local  Cartesian  coordinate  grid  centered  on  Wisconsin. 
Imagery  from  the  AVHRR  3.7  and  11.8pm  thermal 
channels  (  3  and  5 )  is  presented  in  Figures  la  and  1b, 
respectively.  These  images  show  widespread  cirrus 
flowing  southeastward  into  Wisconsin,  just  ahead  of 
an  upper-level  ridge  axis.  Cirrus  are  also  found  over 
southeastern  Wisconsin  near  Madison.  Note  the 
seemingly  more  extensive  detection  of  cirrus  in  the 


Figure  1b.  AVHRR  11.8  pm  Image  for  the  0930  UTC 
28  October  1986  FIRE  IFO 
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11.8  |xm  image,  and  the  apparent  lack  of  cirrus  in  the 
corresponding  regions  of  the  3.7  pm  image,  especially 
thin  cirrus.  This  distinctive  difference  between 
observed  mid-wave  and  long-wave  IR  brightness 
temperatures  enables  AIRC  to  detect  and 
subsequently  analyze  thin  cirrus. 


Having  compiled  the  results  of  the  COz  and  AIRC 
analyses,  a  comparison  was  then  performed  of  the 
independently  retrieved  effective  cloud  altitudes  and 
emissivities. 

4.  COMPARISONS 


All  AVHRR  pixels  in  Figure  1  were  first  checked  for 
cirrus  presence  using  the  detection  algorithm: 

Cirrus:  T3  -  T4  >  0.5  K  ,  (la) 

Thin  Cirrus:  T4  -  Ts  >  0.6  K 

AND  T4  <  273  K,  (1b) 

and  “ Obvious ”  Ci,  Cs: 

T4  <  255  K  ,  (1c) 

where  T3,  T4,  and  Ts  are  the  AVHRR  3.7,  10.7,  and 
11.8  pm  brightness  temperatures,  respectively.  This 
algorithm  is  adapted  from  the  AVHRR  cioud  detection 
scheme  presented  by  Saunders  and  Kriebel  (1988).  If 
any  test  is  satisfied,  cirrus  presence  is  established; 
Figure  2  shows  the  results  of  this  cirrus  screening. 
Inspection  of  animated  sequences  of  GOES  IR 
window  channel  images  verified  these  cirrus  locations 
based  on  their  movement,  and  confirmed'  the  lack  of 
cirrus  in  central  Wisconsin.  Once  identified  as  cirrus, 
each  pixel  triplet  of  AVHRR  brightness  temperatures 
(  T3,  T4,  Ts  )  is  further  analyzed  using  the  AIRC  model 
for  emissivity  and  effective  cloud  layer  altitude. 

Figure  3  contains  the  HIRS  grid  for  the  image 
shown  in  Figure  1.  A  total  of  153  HIRS  pixels  within 
the  Wisconsin  area  were  analyzed  by  Wylie  and 
Menzel  for  effective  cirrus  emissivity  NeC02  and 
effective  cirrus  altitude  zC02  using  the  C02  Slicing 
technique.  Coincident  with  the  HIRS  grid  are  139,649 
AVHRR  LAC  pixels,  each  with  an  AIRC  analysis  of 
cirrus  emissivity  sA|RC  (N  is  assumed  1  for  the  1-km 


Figure  2.  Cirrus  Detection  Results  for  the  Image 
Scene  in  Figure  1.  Bright  Pixels  Denote  Cirrus 


For  each  HIRS  FOV,  a  cirrus  cloud  fraction  N  is 

obtained  from  the  AVHRR  data  using  the  cirrus 

detection  algorithm  specified  by  Eq.  (1).  This  fraction 

is  used  to  separate  the  C02  cirrus  emissivity  ec02  from 

the  effective  emissivity.  The  resultant  11.1  pm 

emissivity  eC02  can  be  compared  to  the  average  AIRC 

10.7  pm  emissivity  eA1RC  of  the  cirrus-filled  pixels  over 

the  same  area.  The  important  point  to  note  here  is 

that  the  AIRC  analysis  allows  for  removal  of  the  cloud 

fraction  dependence  of  the  C02  effective  emissivity 

Nf 

^002* 

Figure  4  contains  a  scatter  plot  of  the  C02  effective 
emissivity  NeC02  vs.  the  AIRC  cirrus  cloud  fraction  N. 
Points  above  the  line  Nec02  =  N  are  spurious  because 
e  is  constrained  to  be  <  1.0,  and  thus  Ne  <  N.  The 
large  majority  of  NeC02  values  above  the  Ne  =  N  line 
represent  C02  opaque  cloud  reports  (  Necoz  =  1  )  in 
areas  where  AIRC  found  little  or  no  cirrus  and  where 
the  C02  radiances  detect  no  transmissive  cloud. 
These  opaque  clouds  were  reported  due  to  an 
overestimate  of  Ts(c  ,  causing  Tg  to  be  flagged  as 
sufficiently  colder  than  the  background  and  hence 
indicative  of  opaque  lower  cloud.  Other  spurious  C02 
points  are  likely  due  to  upper  tropospheric  water  vapor 
in  amounts  that  differ  from  those  the  C02  Slicing 
algorithm  accounted  for. 

In  Table  2  the  C02  emissivities  ec02  are  compared 
to  the  average  of  the  AIRC  emissivities  eAIRC  of  the 
cirrus-filled  pixels  within  corresponding  HIRS  FOVs. 


Figure  3.  AVHRR  10.7  pm  Image  for  the  Analysis 
Area  of  Figure  1,  With  the  HIRS  Grid  Overlay 
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Figure  4.  Scatter  Plot  of  CO,  Effective  Emissivity  NeCQ2 
vs.  AIRC  Cirrus  Fraction  N.  62  HIRS  FOVs  Were 
Determined  “Clear”  by  Both  Techniques 


Std. 

Correlation  Coefficient 

Emis- 

Mean 

Dev’n. 

with  Other 

with  Cloud 

sivity 

e  Value 

Fraction  N 

£C02 

0.46 

0.18 

0.90 

-0.14 

EA!RC 

0.44 

0.13 

0.46 

Std. 

Correlation  Coefficient 

Altitude 

Mean 

Dev’n. 

with  Other 

with  Cloud 

(km) 

(km) 

z  Value 

Fraction  N 

ZC02 

7.2 

1.9 

0.94 

0.06 

ZAIRC 

7.2 

1.0 

" 

0.16 

Table  2.  Statistics  of  Retrieved  Cirrus  Properties 


from  this  analysis  by  restricting  cloud  fraction  to  be 
N  >  0.1  (10%),  and  emissivity  0.1  <  e  <  1.  The 
emissivity  constraint  is  invoked  because  noisy  HIRS 
and/or  AVHRR  radiance  measurements  can  be 
confused  with  very  thin  cirrus,  thus  affecting  the  ability 
of  either  C02  Slicing  or  AIRC  to  detect  and  analyze 
such  clouds.  The  correlation  coefficient  is  0.90 
between  the  two  independently  derived  IR-window 
cirrus  emissivities.  Overall,  the  COz  and  AIRC  models 
calculate  average  cirrus  emissivities  of  0.46  and  0.44, 
respectively. 

Also  listed  in  Table  2  are  the  comparisons  between 
the  C02  effective  cirrus  altitudes  and  the 
corresponding  average  AIRC  altitudes.  These  values 
agree  well,  with  a  correlation  of  0.94.  In  addition,  the 
average  effective  altitude  values  ZC02  and  ZAIRC  for  the 
HIRS  grid  are  equal  at  7.2  km.  The  range  of  retrieved 
cirrus  altitudes  analyzed  by  both  the  C02  and  AIRC 
models  agrees  well  with  coincident  surface-based 
cirrus  lidar  measurements  reported  by  Sassen  et  al. 
(1990). 

Retrieved  emissivity  and  cloud  altitude  were  also 


compared  with  cloud  fraction.  Table  2  shows  the  poor 
correlations  between  e  and  N,  and  between  z  and  N. 
These  indicate  a  lack  of  dependence  of  retrieved 
cirrus  emissivity  on  fraction,  and  a  lack  of  dependence 
of  retrieved  altitude  on  fraction;  such  independence  in 
both  cases  is  physically  reasonable. 

5.  SUMMARY  AND  CLOSING  REMARKS 

Radiative  and  spatial  properties  obtained  by  the 
C02  Slicing  and  AIRC  models  compare  well  for  the 
cirrus  clouds  observed  during  the  0930  UTC  NOAA 
overpass  on  28  October  1986.  The  mean  cirrus 
emissivities  and  altitudes  retrieved  by  the  C02  Slicing 
and  AIRC  models  agree  well  with  each  other;  0.46  : 
0.44  and  7.2  :  7.2  km,  respectively,  with  somewhat 
greater  variance  evident  in  the  C02  values. 

For  other  more  dense  cirrus  cloud  conditions, 
especially  those  of  extensive  precipitating  systems 
where  N  is  large,  C02  Slicing  results  (Wylie  and 
Menzel,  1989)  suggest  that  emissivity  and  fraction 
may  not  be  totally  independent.  Such  cloud  conditions 
will  be  included  in  future  studies.  Also,  the  relationship 
between  cirrus  emissivity  and  altitude  was  investigated 
for  the  28  October  scene  and  no  significant  conclusion 
could  be  reached. 

What  is  indicated  by  the  analysis  of  this 
transmissive  cirrus  scene  is  that  neither  cirrus 
emissivity  nor  altitude  is  dependent  on  cirrus  fraction. 
The  consistent  findings  for  both  the  C02  Slicing  and 
AIRC  models  suggest  that  representative  cirrus  cloud 
properties  are  being  retrieved. 
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